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\q • Abstract 

We study the suppression of high-p T hadron spectra in nuclear collisions, supple- 
menting the perturbative QCD factorized formalism with radiative parton energy loss. 
We find that the nuclear modification factor Raa(pt), which quantifies the degree of 
suppression, is almost pr-independent both for RHIC (in agreement with data) and for 
the LHC This is a consequence of the shape of the partonic p^-spectrum in elementary 
collisions which implies that for the same value of Raa at higher p?, an increasingly 
smaller fraction of parton energy loss is needed. When the values of the time-averaged 
transport coefficient q exceed 5 GeV 2 /fm, Raa(pt) gradually loses its sensitivity to the 
corresponding produced energy density. This is due to particle production in the outer 
corona of the medium, which remains almost unsuppressed even for extreme densities. 
Thus, even for the highest experimentally accessible transverse momentum at the LHC 
and in contrast to jets, the measurement of leading partons via leading hadrons is not 
a penetrating probe of the dense matter, but a rather fragile probe which fragments 
for the opacities reached below the skin of the medium. Relating the transport co- 
efficient to the energy density produced in the collision region, we find and discuss a 
phenomenon reminiscent of the opacity problem of elliptic flow: namely, the interaction 
of the hard parent parton with the medium appears to be much stronger than that 
expected for perturbative interactions of the hard parton with an ideal quark gluon 
plasma. 
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1 Introduction 

In Au+Au collisions at ^/snn = 130 and 200 GeV, experiments at the Relativistic 
Heavy Ion Collider, RHIC, observe a strong, a factor of 4-5, suppression of high-pT 
leading hadron spectra [1, 2, 3, 4, 5, 6]. The centrality dependence of this effect, the 
corresponding suppression in leading back-to-back correlations [7], and the moderate 
enhancement and opposite centrality dependence observed in the same p^-range in 
d+Au collisions [8, 9] are strong indications that the suppression is predominantly a 
final state effect in which the energy of high-p^ particles is degraded depending on 
the in-medium pathlength and the density of the produced medium. Formation time 
arguments [10] imply that for px > 5 — 7 GeV, this energy loss is of partonic origin 
while hadronization occurs outside the partonic medium. This is also supported by the 
disappearance of the particle species dependence of the suppression for pt > 5 — 7 GeV. 
In this high-pr regime, the suppression of hadrons is widely regarded as a "hard probe" 
for the characterization of the energy density of the medium produced in nucleus- 
nucleus collisions. Recent studies, based on medium-induced radiative parton energy 
loss [11, 12, 13, 14, 15, 16], quote time-dependent energy densities which at time r = 0.2 
fm/c after the collision lie a factor 0(100) above the energy density of cold nuclear 
matter [17, 18]. 

In this paper, we calculate the nuclear modification factor for high-pj< hadronic 
spectra in nuclear collisions, 

d 2 N AA /dp T dy 
RAA{PT,y)- 



(T AA ) c d 2 a NN /dp T dy ■ 

Here, (T AA ) C is the standard nuclear overlap function, calculated as the average in 
the measured centrality class. In the absence of nuclear or medium effects, Raa = 1- 
The starting point of our calculation is the collinearly factorized leading order (LO) 
perturbative QCD formalism which reproduces the shape of the hadronic high-pr spec- 
trum d 2 a NN /dpxdy in proton-proton collisions at RHIC energies [19]. Our discussion 
will be limited to pt > 5 GeV. In the calculation of the charged hadron spectrum 
d 2 N AA /dpxdy in this range, we include both initial and final state effects. Nuclear 
modifications to parton distribution functions are taken from the EKS98 parametriza- 
tion [20] and isospin effects are included as in [19]. As a medium effect, we con- 
sider partonic final state energy loss, which is introduced via the quenching weights 
Pf(AE,L,q) calculated in Ref. [21]. The quenching weights will be discussed in more 
detail in section 2.1. They denote the probability that a hard parton / radiates an 
additional medium-induced amount of energy AE if it propagates through a medium 
of length L. The quenching weight depends on the transport coefficient q which is 
determined by the average transverse momentum squared transferred to the projec- 
tile parton per unit pathlength, q = (gf,)/A. It characterizes the energy density of 
the medium, q oc e 3 / 4 [22]. With this input, the hadron spectrum in nucleus- nucleus 



collisions reads schematically 



d^) k+X = E KU f+X ® ^/(A^> ^ Q) ® J>£3(*, a4) • ( 2 ) 

/ 



Here 



dcr (^7 /+X = E /<m( x i» ^ 2 ) ® fj/A(x 2 , Q 2 ) ® £y-/+k (3) 

zjA; 

and fi/A(%,Q 2 ) are the nuclear parton distribution functions and cr„^./ + fc the pertur- 
batively calculable partonic cross sections. 

The paper is organized as follows. In Section 2 we give details of the formalism used 
to compute hadron spectra with parton energy loss. Results for Au+Au collisions at 
v /snn = 200 GeV are presented in Section 3 and predictions for the LHC are given in 
Section 4. In Section 5, we discuss the nuclear modification factor for Au+Au collisions 
at intermediate energy ^ t /s^ = 62.4 GeV, where data are just being analyzed. In 
Section 6, we discuss in detail how the main model parameter of our analysis, the 
time-averaged transport coefficient q, is related to the energy density produced in the 
medium. The main conclusions are summarized in Section 7. The results for RHIC 
and for the LHC are further explored in Appendix A. 

2 Framework 

2.1 Quenching weights 

Under the assumption that multiple gluons are emitted independent of each other, the 
probability that a parton loses an amount AE of its energy by medium-induced gluon 
radiation is 



(4) 
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This probability distribution has been calculated and tabulated [21] for the medium- 
induced gluon energy distribution ujjj^ of a hard parton. [In the present study, we use 
a probability distribution calculated with 0:5 = 0.5 instead of as = 1/3. Since AE oc 
a s q, the dependence on the value of the coupling constant can be largely absorbed in 
a rescaling of q.] The calculation of Ref. [21] was done in two approximations which 
differ in viewing the medium as a source of many soft [12, 13, 14] or a few hard [14, 15] 
momentum transfers, respectively. For the purpose of the present study, the small 
numerical differences between both approximations are not important [21]. We use the 
soft multiple scattering approximation, where the quenching weights P(AE) depend 
on the parton species (light quark or gluon), the in- medium pathlength L and the 



transport coefficient q. For a static medium, q does not show a time-dependence. 
Extension and density of the medium can be expressed in terms of the characteristic 
gluon frequency u c = hqL 2 and the dimensionless quantity R = uj c L. 
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Figure 1 : Continuous part of the quenching weights for light quarks (upper row) and gluons 
(lower row) computed in the multiple soft scattering approximation with as = 0.5. 



In general, the density of the medium decreases strongly in time due to longitudinal 
and transverse expansion. Hence, the transport coefficient should be time-dependent. 
However, the gluon energy distribution uj^ for a time- dependent transport coefficient 
q T agrees with the energy distribution obtained for an equivalent static linear line- 
averaged transport coefficient [24, 25] 



Q = J2 \ rq T (r)dT. 



(5) 



Here, To is the time at which the parton is produced and for r <C L, the To-dependence 
of q is numerically negligible. This allows us to base our analysis on an equivalent static 
scenario and to translate afterwards with Eq. (5) the values for the time-averaged q 
into time- dependent transport coefficients . 

The probability distribution P(AE/u c , R) consists of a discrete and a continuous 
part (see Figs. 1 and 2) 



P(AE/u c , R) = p (R)8(AE/u> c ) + p(AE/cu c , R) . 



(6) 




Figure 2: Discrete part of the quenching weights for light quarks and gluons computed in 
the multiple soft scattering approximation with as = 0.5. 



The continuous part denotes the probability to lose an additional energy AE due to 
medium induced gluon radiation. The discrete contribution p is the finite probability 
for the projectile to escape the collision region of finite size and finite opacity without 
further interaction [24]. The quenching weight P(AE/uj c ,R) is normalized to unity, 



d(AE/u c )P(AE/u c , R) = p (R) 



d(AE/u c )p(AE/u; c , R) 



(7) 



P(AE/u c , R) is a generalized probability. It can take negative values since it describes 
a medium- induced redistribution of gluon radiation which - compared to the gluon 
radiation in the vacuum - can be reduced in some range of phase space [26]. 

The quenching weights are computed in the eikonal limit of very large initial parton 
energy E{. For the case of realistic initial parton energies, finite energy corrections 
have to be taken into account. In the absence of a formalism which includes them a 
priori, this is done by ensuring with kinematic cuts that the probability distribution 
P(AE/u c ,R) has weight only in the range AE < E t [23, 27, 28, 29]. To illustrate 
the theoretical uncertainty associated with finite energy corrections, we shall compare 
results without such corrections to results obtained from the reweighted ('rw') proba- 
bility [28] 



P rw (A£/cu c , R) 



P(AE/cu c ,R) 



J^d(AE/cu c )P(AE/u c ,R) 



9(1 - AE/Ei) . 



We consider the non-reweighted result for the ratio Raa as a lower limit. 



Numerically, we have checked that the effect of reweighting (8) is comparable to the 
finite energy corrections in other approaches [27] where kinematic cuts are imposed 
on the level of the gluon energy distribution by multiplying with the vacuum splitting 
function. 

2.2 Hadron spectra with parton energy loss 

We first discuss the nuclear geometry entering the calculation of the nuclear modifica- 
tion factor (1). The centrality dependence is taken into account by defining an effective 
nucleus A c r for which the number of participants in central A e g + A c q collision equals 
(-/Vpart)c- Here, the number of participants is obtained from an optical Glauber analysis 
with realistic Wood-Saxon (WS) profile for the centrality class c [30]. As far as energy 
loss is concerned, the A c g nucleus is assumed to be cylindrical and uniform, but we 
correct the final spectra for the difference between T A A (0) and T A NS A (0). 

The partons produced at midrapidity y = propagate in the transverse plane. 
Consider a parton which is produced at a position s = (s, <ps) with respect to the 
transverse coordinates and which propagates in the transverse direction given by the 
angle 0^. Its angle of propagation with respect to s is 4>ls — <t>L — <fis- The geometrical 
transverse path length L geom for this parton is 



£gcom(s, (pLs) = S COS (f) LS + yj S 2 COS L5 2 + R Acff - S 2 , (9) 

where Ra cS is the radius of A c q. We assume that significant energy loss occurs only 
in the dense QGP-phase whose estimated lifetime tqgp may be shorter than L ge om- 
Therefore the in-medium path length which enters the quenching weight (5) is chosen 
to be L = min(L geom , L cut ~ tqgp). Varying L cnt will turn out to be a useful tool to 
establish to what extent particle emission is surface dominated. 

We compute the inclusive production of charged particles at midrapidity in the 
framework of factorized leading-order partonic cross sections, using nuclear parton dis- 
tributions, fragmentation functions and quenching weights but no intrinsic transverse 
momenta. As in [19] the renormalization and factorization scales are set to \x — Q — qr 
and fragmentation scale to fiF = Pt, where qr and pt are the transverse momentum 
of the parton and hadron, respectively. We use CTEQ5 parton distributions [31] with 
nuclear effects for A e s from EKS98 [20], and KKP fragmentation functions [32]. The 
p + p reference spectra are obtained as in [19]. 

The inclusive cross section for the production of a parton of flavour / and an initial 
transverse momentum qr i and rapidity yi is given e.g. in Eqs. (4) and (5) of [19]. 
When traversing a medium, the parton will lose a fraction e of its initial energy Ei, 
but it will not change significantly the direction of its trajectory, fa = 4>ls- Therefore, 



we can write for the spectrum of the outgoing partons 

dN AA-+f+x ,1 . _ _ /• , o , da AA ^f +x 



[T A (s)] 2 J deP(e,L,q) fd 2 q Tt dy t 



d 2 sdq Tf d(j)fdyf Jo J " d 2 q Ti dyi 

X-KVf ~ VijK^f - <j>i)8(0r f - (1 - e)q Tl ) 
x0(O <(j)i < 2vr) 

Xe (*^^2^-) e ( fclSarc0S '€)' <10) 

where Ta(s) = , -m \ ®(Ra ~ s2 )' 9o is the lower limit for the initial partonic transverse 
momentum and P(e, L, q) is the quenching weight. In our notation AE = eEi = 
eq Ti coshyi, so 

P(e,L,q) = ^P(AE/uJ c ,R) (11) 

^P(AE/uj c ,R) 
^w(e, L, q) = —. ^ -— - . (12) 

f Q E ^d(AE/cu c )P(AE/u c ,R) 

From Eq. (10) we find that 

max(0, l-qr/q ) <e< 1 *-= (13) 

V s 

and 

dN AA^f+X rR A /-27T /• 1 rfo .A4->/+X 

^PST = 7. dss l ^y* F(£ ' L '' )[T ' 4WI (i-^ dA-fe 



i T .wyi 



X9 (° £ *> £ 2^) 9 l 12 " 1 £ arC ° Sh 



€' ■ <"> 



with q Ti = q Tf /(l - e) and y» = y/. 

We define z as the fraction of the energy of the parent parton carried by the leading 
hadron, z = E^/Ef. For a hadron with transverse momentum p? we obtain [19] 

, N AA^h+X r j dN AA-^f+X 

Ki^.J^y^j^D^^JL^-, (15) 

qT f = — J(rtvr,y), Uf = arsinh I sinhy I (16) 

(17) 



dp^dy 
where 

and 



z 




\Pt 


J(m T ,y) = 


(- 


2 \ -V2 


m\ cosh 2 y J 



The integration region for z is 

2m,T 

7^ 



cosh y < z < 1 , 



since we do not attempt to explore the region p? < qo- The cms-energy dependent K- 
factors applicable within our framework were determined in [19]. They are important 
for the overall normalization of the spectrum but they cancel in the nuclear modification 
factor Raa discussed in the next sections. 

Our formalism is equivalent to defining medium-modified fragmentation function 
[24, 27, 33] 

*>£?(*, *40 = J deP(e)^L_D f ^ h (^,^) • (19) 

Neglecting the masses, Eq. (15) simplifies [23] for the case of small AE and z = 
E h /E f ~ 1 to 



i"nj 



dprdy J dprdy 

Here, dN^ va ~ c ^ / dprdy is the spectrum of hadrons in the case of no medium (i.e. it 
coincides with the proton-proton spectrum in the absence of initial state effects as 
shadowing etc.). The suppression computed with (20) gives a rather good approxima- 
tion to the one computed with the full formula (15). 

3 Results for RHIC at y^ = 200 GeV 

Based on the above formalism, we have studied the nuclear modification factor Raa{pt) 
for charged hadrons, h = (h + + hr)j1 in Au+Au collisions at ^Js^n = 200 GeV. 
We include charged pions, kaons and (anti)protons in the sum of charged particles. 
Our discussion is limited to the — 5% most central collisions. This centrality class 
corresponds to the head-on collision of an effective nucleus of size A c g = 181. 

In Fig. 3, we compare data from RHIC with the nuclear modification factor calculated 
for different values of the transport coefficient q. This figure is for a medium which 
exists for at most r ~ L cut = 5 fm in the dense partonic phase. We find that even in the 
absence of final state parton energy loss, q = 0, the nuclear modification factor shows 
a reduction from the baseline Raa = 1 for p? > 13 GeV. The reason is that in this pt 
region, the incoming parton momentum fractions x ~ Pl "4 sampled by the hadronic 
spectrum lie in the so-called EMC region where nuclear parton distribution functions 
decrease faster than free proton parton distribution functions with increasing x [20]. 
This decrease of the incoming parton flux translates into a pr-dependent decrease of 
the nuclear modification factor Raa- As the density of the medium (or equivalently the 



transport coefficient) is increased, the nuclear modification factor is reduced further. 
This is caused by parton energy loss in the final state. Interestingly, the resulting 
nuclear modification factor Raa is almost pr-independent and its dependence on q 
becomes weaker as q increases. We now discuss these two features in more detail: 
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Figure 3: (Colour online) The nuclear modification factor Raa for charged hadrons h = 
(h + + h~)/2 in the — 5% most central Au+Au collisions at ^/snn = 200 GeV. Different 
lines are for a maximal lifetime of the medium r ~ L cu t=5 fm and different values of the 
transport coefficient q. The shaded region between the curves calculated with reweighted 
(solid curves) and non-reweighted (dotted-dashed curves) quenching weights is indicative of 
uncertainties related to finite energy corrections. For comparison, the most central Raa data 
from the RHIC collaborations are shown (with statistical errors only): — 5% STAR data 
for |t/| < 0.5 [4], - 10% PHENIX data for |r/| < 0.35 [2], - 10% BRAHMS data for 77 = 
[6] and - 6% PHOBOS data for 0.2 < rj < 1.4 [5]. 



To understand the approximate pr-independence of Raa, we study the sensitivity 
of Raa to the pr-dependent shape of the partonic cross section entering (2). (For 
qualitative results we assume z ~ 1, which is a fairly good approximation for quark 
dominated hadron production at RHIC in the large-p^ region of our interest.) We 
consider first a partonic cross section which has a powerlaw tail da /dp? ~ l/pj 1 , n 
constant. If n is sufficiently large, then the medium-modified spectrum can be obtained 
by shifting the vacuum spectrum by S(jpt) ~ \/pT/ n [23], and one finds Raa(pt) ~ 
(l + c/y/npr) n where c = y/2'KoP' l xo c . In this case, Raa approaches unity in the limit 



of large px, except for possible modifications in the incoming parton distributions. In 
contrast, if the partonic spectrum is exponential, one finds a pr-independent nuclear 
modification factor 



Raa(pi 



d(AE)P(AE) 



exp{-q(pr + AE)} 
exp{— clpt} 



d{AE)P(AE) exp{-aAE} . 

(21) 

For the spectra of produced quarks and gluons computed from Eq. (15) at RHIC ener- 
gies, one finds that the logarithmic slope n of the partonic p^-spectrum increases with 
increasing px- This makes parton energy loss more effective in reducing Raa at high 
Pt- Moreover, due to kinematic boundary effects, the underlying partonic spectrum at 
RHIC can be parametrized by an exponential exp{— apx} for px > 30 GeV. As a con- 
sequence of (21), one thus finds - in agreement with Fig. 3 - an almost p^-independent 
Raa at large pt > 30 GeV. [Although this is not within experimental reach, we note 
that close to the kinematical limit px ~ 100 GeV, the slope decreases even faster.] The 
remaining weak p^-dependence can be attributed to the parton species dependence of 
the quenching weights and to the nuclear dependences of the incoming parton distribu- 
tion function which turn out to further reduce Raa with increasing px- In conclusion, 
the onset of a perturbative px powerlaw (at px > 5 GeV say), and its return to an 
exponential spectrum at px > 30 GeV leave a finite window at RHIC energies. Within 
this window, the gradually increasing slope of the partonic spectrum implies an almost 
p-r-independent nuclear modification factor. 
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Figure 4: (Colour online) Raa for charged hadrons h = (h + + h )/2 in the - 
central Au+Au collisions at y/SNN = 200 GeV as a function of q for p^=10 GeV. 
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We now turn to the dependence of Raa on q. Partons which traverse a larger in- 
medium pathlength tend to lose more energy. As a consequence, the dominant fraction 
of leading partons comes from the outer skin of the two-dimensional transverse overlap 
of the colliding nuclei. We refer to the outer region of the medium as "skin" rather than 
"surface" since it has a finite three-volume. In Appendix A, we give further support 
for this dominance of skin-emission for realistic densities at RHIC and at the LHC. 
The dominance of skin-emission appears to limit the sensitivity of Raa on the density 
of the medium. Indeed, we find that for sufficiently large but realistic densities of the 
medium, the thickness of the outer skin changes only weakly with increasing density, 
see Fig. 4. Thus the nuclear modification factor is a sensitive measure of the density 
of the medium up to q < 4 GeV 2 /fm, but it loses this sensitivity for higher values of q. 

From Fig. 3 and 4 we conclude that the suppression of high-pT hadrons measured 
in v /snn = 200 GeV Au+Au collisions at RHIC requires a time-averaged transport 
coefficient q > 5 GeV 2 /fm and favours q = 5 — 15 GeV 2 /fm. 

4 Extrapolation to the LHC 

Within the formalism of section 2, the choice of the transport coefficient q is the main 
unknown for extrapolating our results from RHIC to Pb+Pb collisions at the LHC 
energy ^/snn = 5500 GeV. One expects that this transport coefficient grows linearly 
with the number density, rij, of the medium [22, 25], and hence is proportional to the 
multiplicity per unit rapidity 

dN ch 
?oc^-. (22) 

dy 

For the LHC, model predictions of the event multiplicity vary by a factor ~ 4 at 
least [34], but the corresponding uncertainty for Raa m ay be significantly smaller 
since Raa changes only weakly for large values of q. The second unknown in the 
calculation of Raa is the maximal lifetime of the dense medium, denoted by L cnt . 
We expect, however, that this choice is numerically unimportant since L cnt at LHC is 
typically larger than the thickness of the effective skin of the medium through which 
hard partons can penetrate without significant energy loss. 

To be specific, we use for the LHC L cut = tq GP ~ 10 fm consistent with hydrody- 
namical modeling of the time evolution of the medium [35]. We find from an optical 
Glauber calculation that the — 5% most central Pb+Pb events correspond to the 
head-on collision of nuclei with A e s = 193. The density estimate is obtained from the 
EKRT parametrization [36] for which rii ~ A 0383 ( v /i) - 574 , consistent with the growth 
of the event multiplicity measured so far. With this input, the transport coefficients 
at RHIC and the LHC are related like ^lhc — 6.8gRHio F° r a typical RHIC value, 
<7rhic — 10 GeV 2 /fm, we find that ^lhc — 68 GeV 2 /fm which is one of the parameters 
used in Fig. 5. In models which predict a smaller multiplicity increase from RHIC to 



10 



LHC [34], the suppression is expected to lie in between the bands shown in Fig. 5 for 
g RHIC = 10 GeV 2 /fm and q LHC = 68 GeV 2 /fm. 



Qc; 



1.2 

1.1 
1.0 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



1 1 1 1 1 1 1 1 1 1 1 1 


l 


- f — — 


q [GeV7fm]= 


- 


Ojno medium 


- 


1 , rw 

— " ! ~ 


- ^rf^""-"" — 


'_ /^ vVw = 5500 GeV, {h + +h)l2 


- 


A e£f =193(0-5%) 


- 


~i L cm =10.0fm 


10, rw ~ 

10 


^-^^^ 


68, rw 


/ 


68 




— 


[,i,i,i,i,i, 


1,1,1, 



50 



100 150 200 250 300 350 400 450 500 



Figure 5: (Colour online) The nuclear modification factor for charged hadrons h = {h + + 
h~)/2 in — 5% most central Pb+Pb collisions at the LHC, calculated for L cut =10 fm and 
various values of the transport coefficient. The largest value <Jlhc = 68 GeV /fm lies in the 
range favoured by current multiplicity estimates. See the text for more details. 

Based on estimates of the luminosity in Pb+Pb collisions at the LHC and the 
hadronic cross section, one expects that the nuclear modification factor will be mea- 
sured up to pr ~ 100 — 150 GeV. Not only in this p^-range but even above, the nuclear 
modification factor turns out to be almost pr-independent for <?lhc = 68 GeV 2 /fm. 
Moreover, it is numerically comparable to Raa at RHIC, calculated for a much smaller 
transport coefficient. This can be understood as follows: the contribution from the 
discrete part of the quenching weight to Raa decreases with the increasing transport 
coefficient from RHIC to the LHC (see Fig. 2 and Appendix A). At the same time, 
however, since the partonic transverse momentum spectra are steeper at RHIC than 
at the LHC, the suppression of Raa induced by the continuous part of the quenching 
weight (parton energy loss) is still larger at RHIC than at the LHC (see discussion in 
Appendix A). If the density of the medium created at the LHC is significantly reduced 
(<7lhc — 10 GeV 2 /fm), then one finds a weak but experimentally accessible increase 
of Raa(pt) with pt- The reason is that for sufficiently low density, the probability 
that partons from deeper layers inside the medium escape with small energy loss grows 
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significantly with increasing pt- This effect is even more pronounced for very small 
densities of order ^lhc = 1 GeV 2 /fm. 

From this exercise, we expect that up to the highest transverse momenta accessible 
at the LHC, leading hadron spectra will show a strong suppression in nuclear collisions. 
This suppression is likely to result in an almost pr-hrdependent nuclear modification 
factor. For the reasons given above and since particle emission from the outer skin 
of the medium remains essentially unsuppressed even for extreme densities, we expect 
that the numerical value for Raa at the LHC is comparable to the one at RHIC. 

5 The nuclear modification factor at y / SNN = 62.4 
GeV 

Experiments at the Relativistic Heavy Ion Collider RHIC have just completed data 
taking for Au+Au collisions at an intermediate center of mass energy of ^/snn — 62.4 
GeV. In Fig. 6, we give results for the corresponding nuclear modification factor in 
the present formalism. The main model input is again the value of the time-averaged 
transport coefficient q defined in (5). Based on the multiplicity scaling of the trans- 
port coefficient in (22) and the value q ~ 5 — 10 GeV 2 /fm favoured for RHIC data at 
^/snn = 200 GeV, we chose q ~ 3 — 5 GeV 2 /fm at ^/snn = 62.4 GeV. For the second 
model-dependent input, the life-time L cut of the medium, we chose a slightly reduced 
value L cnt = 4 fm. However the dependence of our results on L cnt turns out to be neg- 
ligible for sufficiently large L cnt - For medium densities favoured by multiplicity scaling 
the suppression is pronounced. The results shown in Fig. 6 are, within theoretical un- 
certainties, consistent with the results published in several recent studies [37, 38, 39]. 

Some caveats should be kept in mind. First, at intermediate pt < 5 — 7 GeV, other 
physical mechanisms can provide sizeable contributions to Raa- In particular, the 
anomalous baryon enhancement at intermediate pt [40, 41] and the Cronin effect [42] 
are observed to increase Raa- Both contributions may be more pronounced on a 
spectrum with steeper slope, i.e. for smaller ^/snn- Second, the smaller the hadronic 
Pt-, the larger the theoretical uncertainty in calculating parton energy loss. Moreover, 
below pt ~ 4 GeV, the hadron spectra at ^/snn = 62.4 GeV deviate significantly from 
the shape calculated by the collinear factorized pQCD formalism [19]. 

While completing this work, the first experimental data on the nuclear modification 
factor in Au+Au collisions at ^/snn = 62.4 GeV appeared [43]. These data are for 
0.2 < r] < 1.4 and p T < 4 GeV with R A a(pt ~ 4GeV) ~ 0.8 ± 0.1 for the - 6% most 
central bin. Other preliminary data [44] indicate a significant further drop of Raa for 
4 < pt < 6.5 GeV. For the highest pt, these data are close to the calculation shown in 
Fig. 6. 
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Figure 6: (Colour online) The nuclear modification factor for charged hadrons h = {h + + 
h~)/2 in — 5% most central Au+Au collisions at y/s = 62.4 GeV. See the text for details. 

6 Relating the Transport Coefficient to Energy Den- 
sity 

In section 3, we have shown that RHIC data favour a time-averaged transport coefficient 
q ~ 10 GeV /fm. This parameter as well as the time-dependent transport coefficient 
entering eq. (5) scales with the energy density of the medium like [22] 



QAr) 



ce^'Hr) 



(23) 



Here, c is a proportionality constant which can be calculated for specific models of the 
medium. In particular, for the model of an ideal quark gluon plasma whose constituents 
interact perturbatively with the hard parton, we extract this proportionality constant 
from Figure 3 of Ref. [22] 

cgg*«2. (24) 

Alternatively, one can determine the constant c from experiment if one has independent 
information about transport coefficient and energy density. From the measured time- 
averaged transport coefficient q, we find via the dynamical scaling law (5) the transport 
coefficient for an initial time tq, 



Qt{To) 



n I L 
To 



(25) 
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This expression is for an expanding medium with q T (r) = q T (T"o) (tq/ L) a . We have 
assumed To <^ L. The expansion parameter a is unity for a one-dimensional Bjorken 
expansion and is expected to stay close to unity for realistic expansion scenarios. From 
this we find 

9t(t ) q 2-a f L\ 

(26) 



e(r )3/4 e ( ro )3/4 2 \r 

For a typical initial time To ~ 0.2 fm/c, the average in-medium pathlength can be 
expected to be ~ 10 r = 2 fm (see Appendix A). The energy density averaged over a 
uniform transverse profile at initial time To can be as large as e(r ) ~ (9(100 j^f), as 
suggested by the pQCD+saturation+hydrodynamics model [36] which correctly pre- 
dicted the final multiplicities and transverse energies in central collisions at RHIC 
[30]. This implies that realistic scenarios lie in the parameter range e(ro) < 100 ^x, 
L ~ 10 ro, 0.75 < a < 1.5, for which we find 

c> 8 . . . 19 . (27) 

Thus, for the entire realistic parameter range, we find a proportionality constant which 
is a factor ~ 4 — 5 larger than the perturbative estimate. We note that q is the product 
of the density of scattering centers times their elastic scattering cross section [22] . Thus, 
the factor ~ 4 — 5 discrepancy observed here is reminiscent of the opacity problem in 
elliptic flow, where also scattering cross sections much larger than the perturbative 
ones are favoured by the data [45]. 

7 Conclusions 

In this paper we have studied within the LO pQCD formalism [19] the effect of medium- 
induced parton energy loss on charged hadron spectra. At high pt we find generically 
an almost p-r-independent nuclear modification factor at both LHC and RHIC energies. 
This is a direct consequence of the interplay of medium-induced energy loss and of a 
partonic p-r-spectrum whose slope increases with increasing p?. In this approach, an 
almost pT-independent Raa is obtained without fine-tuning of nuclear effects (such as 
Cronin or anomalous baryon enhancement). 

For central Au+Au collisions at ^s^n = 200 GeV, we find a satisfactory agreement 
with experiment for a transport coefficient in the range q ~ 5 — 15 GeV 2 /fm. For 
these large values of q, the nuclear modification factor is only weakly dependent on 
the transport coefficient, which makes it difficult to constrain the energy density with 
-RauAu, see Fig. 4. On the other hand, this weakened sensitivity of Raa makes our 
prediction for the nuclear suppression at the LHC more stable. Irrespective of the 
uncertainty in the factor by which the density of the medium increases from RHIC to 
the LHC, the nuclear modification factor is expected to lie in a narrow band -RpbPb ~ 

14 



0.2±0.1. This and the almost pr-independent slope of -RpbPb are the main predictions of 
the present work for Pb+Pb collisions at the LHC. Thus, our calculations indicate that 
even in a kinematic regime pt ~ 100 GeV where the same energy loss formalism results 
in jet cross sections which are almost unaffected by the medium [28], the suppression 
of leading hadron spectra is still sizeable. This points to the importance of resolving 
the internal structure of jets, i.e. the jet shape and jet multiplicity distributions, at 
the LHC in order to access the medium-dependence. 

In section 6, we have inverted the standard logic of the analysis of nuclear modi- 
fication factors. Rather than to determine the energy density e from the measured 
transport coefficient, we extract the proportionality factor c = qe~ 3 ' 4 for realistic up- 
per limit estimates of the energy density. This proportionality factor characterizes the 
strength of the interaction between the hard (test) parton and the medium. It turns 
out to be at least a factor « 4 — 5 larger than expected from perturbative estimates 
[see discussion of eq. (24)]. Although this perturbative estimate (24) is itself subject 
to significant uncertainties, we do not expect that this is the only source of such a large 
discrepancy. Here, we note that our finding is reminiscent of parton cascade studies of 
elliptic flow [45] where perturbative partonic cross sections have to be scaled by large 
factors in order to match the size of the collective evolution. Both effects may have 
a common origin since both concern the interaction strength of partonic degrees of 
freedom with the medium produced in the collision. 

Note added: While completing this work, a calculation of nuclear modification 
factors in a closely related model appeared [46] . This study supports the approximate 
^-independence of the nuclear modification factor at RHIC and at the LHC. Moreover, 
this model calculation reproduces the observed centrality dependence of Raa an d of 
the suppression of the away-side jet-like correlations. That the centrality dependence 
of suppressed high-p^ hadroproduction follows naturally, once the suppression for the 
most central bin is fixed, was also found in Ref. [47]. 
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A Further studies of the dependences of Raa 

We have argued in section 3 that for the densities attained at RHIC, leading hadropro- 
duction is limited to the outer skin of the collision region. Here, we further support 
this argument by calculating the dependence of Raa on L cut . In particular, Fig. 7 
establishes that a reduction of the maximum in-medium pathlength from 6 fm to 4 fm 
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leaves the nuclear modification factor almost unaffected. The slight sensitivity of Raa 
to L cut =3 fm indicates that a small contribution comes from partons which traverse 
3 — 4 fm of matter with transport coefficient q = 10 GeV 2 /fm. Thus, essentially all 
partons with L > 3 fm, as well as a fraction of the partons with L < 3 fm do not 
contribute to the leading hadron spectrum. The average geometrical path length is 
larger than L = 3 fm, it is (L geom ) = -§^Ra = 5.2 fm for a cylindrical uniform nucleus 
(we use A c g = 181 and Ra cS = 6.2). Thus, much less than half of the partons will 
succeed in escaping the medium with sufficient energy, in agreement with the observed 
value for -RauAu at RHIC. 
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Figure 7: (Colour online) The nuclear modification factor Raa for charged hadrons h = (h + + 
h~)/2 in the — 5% most central Au+Au collisions at ^/snn = 200 GeV. The solid curves 
(dotted-dashed curves) are the reweighted (non-reweighted) results with q = 10 GeV /fm 
and L cut =3, 4, 5 and 6 fm. The data are the same as in Fig. 3. 

Another illustration of skin-dominated leading hadroproduction in nucleus-nucleus 
collision is to compare Raa calculated for several fixed in- medium pathlengths L = 1...6 
fm, see Fig. 8. For a small transport coefficient, q = 1 GeV 2 /fm, Raa decreases 
significantly even for L=4, 5, 6 fm, indicating that leading hadrons can originate from 
partons with significant in-medium pathlength. Fig. 8 is consistent with the statement 
in [21] that the suppression of hadrons at RHIC can be accounted for by q = 1 GeV 2 /fm 
and L ~ 6 fm. In the present study, we find a significantly larger transport coefficient 
since the relevant in-medium pathlength turns out to be much smaller. Consistent with 
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the statement made above, if the transport coefficient is increased to q = 10 GeV 2 /fm, 
partons with in-medium pathlength more than 3 fm make a negligible contribution to 
the nuclear modification factor, see Fig. 8. This is consistent with the model discussed 
in Ref. [48]. Effectively, with q = 10 GeV 2 /fm the results from RHIC can be reproduced 
with a fixed path length L ~ 2 fm, which is clearly smaller than {L geom ) = 5.2 fm (see 
also [29]). For L ~ 2 fm, the average momentum transfer acquired by the hard partons 
is of the order y/qL ~ 4.5 GeV. Thus, the direction of partons with q T . < 5 GeV 
is randomized; they become part of the medium. This suggests that hadrons with 
Pt ^ 2 — 3 GeV, which form the bulk of the multiplicity, originate from decoupling of 
the strongly interacting system. This observation is consistent with the results in [49] . 
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Figure 8: (Colour online) The non-reweighted Raa results at RHIC computed using constant 
in-medium path length L for q = 1 GeV /fm (left panel) and q = 10 GeV /fm (right panel). 
See the text for further details. The data are the same as in Fig. 3. 

We finally check to what extent the effects of parton energy loss on Raa are due to 
the discrete or the continuous part of the quenching weight (6). The discrete part p 
depends for fixed L and q only on the parton species, see Fig. 2. The non-reweighted 
contribution of p to Raa therefore interpolates between the values pf, uon and £>o Uar 
as a function of px- In Fig. 9 we fix L = 3 fm, and show separately the discrete and 
continuous contribution to the nuclear modification factor Raa a t RHIC for q — 1 
and 10 GeV 2 /fm. No other nuclear effects are included in Fig. 9. The high-p^ hadron 
spectrum at RHIC is dominated by quark production, and thus the non-reweighted 
discrete component of Raa is very close to £>g Uar • 

At the LHC, for q — 1 and 10 GeV 2 /fm, the suppression is weaker than at RHIC, 
see Fig. 10 and Sec. 4. The contribution of the discrete part to Raa is smaller than at 
RHIC since a significant part of the hadron production is now of gluonic origin. [The 
slight decrease of the discrete contribution at pt > 200 GeV reflects the relative shape 
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Figure 9: (Colour online) The relative contributions of the discrete and continuous parts of 
the quenching weight (6) to the nuclear modification factor Raa at RHIC for fixed in-medium 
path length L = 3 fm and for q = 1 GeV /fm (left panel) and q = 10 GeV /fm (right panel), 
respectively. The data are the same as in Fig. 3. 



of the KKP quark and gluon fragmentation functions at very high z.\ The dominant 
contribution to Raa comes from the continuous part of the quenching weight. Similar 
to RHIC, the increasing slope of the partonic spectra at high p? tames the growth of 
Raa- 
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Figure 10: Same as in Fig. 9 for the LHC. 
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